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Using the method of locally equilibrium statistical operator we consider the thermalized rela- 
tivistic quantum fields in an oscillatory trap. We compare this thermal picture of the confined boson 
gas with non-relativistic model of independent factorized sources. We find that they are equivalent 
in the limit of very large effective sizes R, more exactly, when the Compton wave length 1/m and 
thermal wave length 1/VmT are much smaller than R . Under this conditions we study the in- 
fluence of Bose condensation in finite volumes on the structure of the Wigner function, momentum 
spectra and correlation function. 

I. INTRODUCTION 

In heavy ion collisions at RHIC and LHC energies quasi-macroscopical systems containing 10 4 — 10 5 particles are 
expected to be created. If phase-space densities of such systems at a pre-decaying stage will be high enough, one can 
observe multi-boson effects enhancing the production of pions with low relative momenta, softening their spectra and 
modifying correlation functions. One can even hope to observe new interesting phenomena like boson condensation 
in certain kinematic regions with a large pion density in the 6-dimensional phase space: / = (27r) 3 <i 6 n/<i 3 p<i 3 x >~ 1 
(see, e.g., §||). 

Generally, the account of the multi-boson effects is extremely difficult task. Even on the neglection of particle 
interactions in the final state the requirement of the BE symmetrization leads to severe numerical problems which 
increase factorially with the number of produced bosons [l],|^]. In such situation, it is important that there exists a 
simple analytically solvable models || allowing for a study of the characteristic features of the multi-boson systems. 
Actually, there are two basic methods that are presently used to describe multi-boson systems in finite (small) volumes 
typical for A+A collisions. First one is maximally closed to the procedure of computer simulations of multi- 

boson effects. We call it: the model of independent factorized sources, or MIFS. It implies the independent emission 
of 'probed' non-identical particles with factorized Wigner functions 

n 

D n (pi,xi-...;p n ,x n ) = Y\_D{pi,Xi). (1) 

i=l 

The basic function in such an approach is the 'probed' Wigner distribution D(p,x) that is chosen usually in the 
thermal-like (Boltzmann) form: 

2 2 

Dip, x) = r !— T ^ exp( — ^— - ■^)6(t). (2) 

yy ' ' (27ri?A) 3 yK 2A 2 2R 2> y ' y ' 

We normalize it to a mean multiplicity r\. In typical cases the input numbers n of particles are described by the 
Poissonian distribution 

P(n) = e~VM (3) 

In order to obtain an output ('true') particle number distribution as well as single and two-particle momentum 
spectra, the special procedure of "switching on" of the Bose-statistics is applied. In other words, the symmetrization 
of the 'probe' n-particle normalized amplitudes A n {pi\ is provided for. Then the final multiplicity distribution is 
P s {n) = P{n)N[Al] , where N[A%} is the normalized weight due to the symmetrization. 
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The other approach [|8|]9| deals with true boson statistics from very beginning and is based on the density matrix 
formalism. In particular, the specific p— matrix ansatz in the wave packet formalism was proposed J^] to get (and 
develop) MIFS results. 

To treat the results of these approaches one can often use the statistical thermodynamics language like "rare gas", 
"Bose condensate", etc. At the same time no systematic consideration based on the thermal matrix density for such a 
kind of system in finite volumes has been done. For thermal relativistic essentially finite (small) systems the problem is 
mathematically rather complicated, however for large enough ones it can be solved analytically. We will demonstrate 
it here using the method of statistical operator. 



II. DENSITY MATRIX FOR LOCALLY EQUILIBRIUM SYSTEMS 



According to the definition the statistical operator is 

p=e~ S . (4) 

For locally equilibrium systems the entropy S has to be maximized under additional conditions fixing the average 
densities of energy e(x), momentum p(x), charge q(x). Systems are considered on some hypersurface a : da v = dan v 
with a time-like normal vector n". In the relativistic covariant form the conditions look like 

pP(x) = (e(x), p(x)) = (n u (x)f^(x)) , q(x) = (n u (x)J^x)) (5) 
where (...) means the average with the statistical operator p in (Q). The result for entropy operator is then [jn^, |ll[] 
S = S(a) = $(tr) + / da v (x)^(x)f^(x) - f d* v (x)n{x) J v (x), (6) 



where $(cr) = In Spexp{J da n u {x)(3 il (x)T tLU (x)} is Masier-Planck functional, (3^(x) and p,(x) are Lagrange multipli- 
ers. For real one-component free scalar field we will use the current of particle number density J v {x) = ^ + \x) d v 
ip^(x) where and <p(' are the positive and negative field components. The energy-momentum tensor has the 
standard form. 

Then for the system which has no internal flows and is in local equilibrium state on the hypersurface t = 0: 



n" = (1, 0), at the same temperature T : /3^(x) = (^,0) the density matrix looks like: 
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9 = ~z cxp 



-0 J d 3 pp°ala p + ^-y 3 J d 3 xp(x)^^(k + k' Q )e-^-^ala kl 
To guarantee the maximal closeness to the MIFS we will use the chemical potential in the form: 



(7) 



The description of inclusive spectra and correlations for a multiparticle production is based on a computation of 
the averages 

= ^ ^ ' P ° lP ° 2 dp!dp? = ( a pi a i a Pi a p 2 ) i etc >- (9) 
For Gaussian-type operator like (]?]) the thermal Wick theorem takes place: 

(a+ a+Opiflpa) = (a^a Pl ) (a+ 2 a P2 ) + (a+a P2 ) (a+ a Pl ) . (10) 
To find the averages {a pi a P2 ^ the Gaudin's method [Q is used [^3|. As a result we have: 

OC r. 

(a+a P2 ) G "bi'P2); G n (pi,p 2 ) = J d 3 kG n - 1 (p 1 ,k)G 1 {k,p 2 ). (11) 



71=1 



The basic function G\(p\,p2) will be defined below. For this aim, first the commutation relation with entropy operator 
have to be considered: 
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[a(p),S] = J d 3 kM(p,k)a(k) M(p,k) « M<°>(p,fc) + O(^) + °(^) ( 12 ) 
where, using (|7|) and (g), we have 

MW(pi,p2) = 0p° 2 S( Pl - pa) - | (I^e-fPr^X,). (13) 



The basic function Gi(pi,p 2 ) is then defined as follows: 

Gl(pi, P2 ) = Y / ^-r M n(PuP2), (14) 
n=0 

where 

M)(Pi, P2) = S(pi - p 2 ), Mi (pi, p 2 ) = M (pi, p 2 ); 

(15) 

Afn(Pi,Pa)= / d 3 kM n -i(pi,k)M 1 (k ) p 2 ). 



III. THE COMPARISON WITH MIFS. 

At m 2 R 2 ,mf3- 1 R 2 > 1 we find, neglecting the terms 1/pq-R 4 and P 2 /Po r4 in Eqs.(|l^), 
Gi(p llP2 ) = G£ 0) (pi,pa) + G^(px,p 2 ) + 0(-^, ^) « 

^/Ae^expf-^"-^)))/ (16) 
Here p = (pi + p2)/2; q = pi — P2- At m 2 R 2 , A 2 R 2 — » oo in the non-relativistic limit the Wigner distribution is 

£ p2 j.2 

D ^ p) ^j2^ eXP ^h-^ (17) 

where A 2 = mT and fugacity £ = e M , /I = /1 — m. Comparing Eqs. ( [T7| ) and (Q) we found that in the non-relativistic 
approach in the limit of very large emission volumes the thermal density matrix (J7]) with the chemical potential (|^) 
and the fugacity £ = r)/(AR ) 3 is equivalent to MIFS. The physical reasons are the following. First, if the wave-length 
of the quanta is much less that the system size, the particle does not 'feel' the finiteness of the systems and the 'probe' 
distribution will be the same as in the thermodynamic limit, i.e., the Boltzmann-like one. Second, for large systems 
the assumption of an independent emission of non-interacting particles is natural: the average distance between any 
two particles is large comparing with its wave-length (or with the size of the wave packet). 

At the same time Eq.(Q) indicates that for essentially small system's sizes compared with quanta wave-lengths, the 
single-particle locally equilibrium distribution D(x,p) has no more simple Boltzmann-like form (^). Some distortion 
terms have relativistic nature and are essential when the system size is close to the Compton wave-length, R ~ l/m. 
The others can take place even if Rm 3> 1, depending on the ratio between size R and thermal wave-length 1/ V mT. 
They are of quantum nature and appear when one describes the thermal equilibrium of the quanta with an average de 
Broglie wave length that is larger than the system size. In general, according to Eq.(|l^) one can expect the reduction 
of soft quanta in thermal model in comparison with the MIFS Boltzmann anzats while the distributions of "hard" 
quanta coincide in both approaches. 
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IV. THE WIGNER FUNCTIONS AND SPECTRA IN MULTI-BOSON SYSTEMS 



We will consider the limiting behavior of the functions G„(pi, P2) ( |ll[) at small n << AR and large n ^> AR at 
RA ^> 1 in the non-relativistic case. Neglecting the terms 1/A 2 R 2 and l/m 2 R 2 in Eq. ( ]T6| ) we can put C?i(pi,P2) = 
Gi (pi,p2) and found the optimal tailing of the two limiting forms at the point nt = RA. Then at RA 1 one can 
express the basic operator average ( pi] ) through the Wigner functions of the Bose gas (g) and the Bose condensate 
(c) (we will see the correspondence later) as follows: 

T «0 =f]G n (p 1 ,p2) « / ^xe^{f g {p,x) +/ c (p,a;)), (18) 

7 2 . 1 J 



Pi 



where 



fg(P,x 



(19) 



and 



(2^) 3 ^ 1 exp(p 2 /2A 2 )-l ' 
£, = £cxp(-x 2 /2i? 2 ) = e^^-sss) exp(-x 2 /2i? 2 ) 

/ c (p,z) = J^«Aj_ exp( _^ p2 _ A a 

1 — £ 7T A it 

We have used here the following transformation: 

Note here that the critical value of the nonrelativistic chemical potential = ^ s positive because gas is not 
ideal (effectively the gas interacts with external field confining it in a finite volume). 

In the thermodynamic limit: volume V — > 00 (RA — > 00) at fixed temperature T and at fixed particle density 
n = N/V the phase space density in the central part is 



where § a (z) — ^ZfcLi w • ^ ne nrs t term corresponds to the ordinary Bose- Einstein gas contribution. Let us introduce 
associated with it at £ = 1 the critical density: n cr — ( 2 J> 3 / 2 $3/2(1)- It is easy to see that in order to move to the 
thermodynamic limit at fixed density n > n cr the parameter £ at large enough RA has to tend to unity as: 

^ = 1 ~ [n-n cr )(-KRAf/ 2 ' ^ 

The single particle inclusive distribution in the thermodynamic limit has the form: 

00 

n (1) (p) = I d 3 xf(x, p) = n g (p) + n c (p) « 

— OO 

j^m* 3 / 2 &-^) + X^^i? 3 (^A)- 3 / 2 exp(-f p 2 ) (24) 
RA ^ jfr I d3 * l^i^i + (^ A ) 3/2 (- ^ "Mp) 



Two terms in the last equality of Eq. ([24]) are associated with the classical Bose gas and the condensate in thermody- 
namic limit at the densities n > n cr . They describe Bose gas and Bose condensate for finite systems at the condition 

RA » 1. Note that for Bose condensate the momentum distribution is much more narrow than for Bose gas, J 
« A. 
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V. THE INTERFEROMETRY OF MULTI-BOSON SOURCES 



The calculation of two-particle inclusive spectra is based on Eqs.(||), jlO|), (pf). One can easily find that 

-i (a+a p2 ) = n g (p) exp(-q 2 i? 2 /2) + n c {p) exp(-q 2 i? 2 /2), (25) 

where 

R 2 g = R 2 $ 5/2 &-^)/<i> 3/2 &-^) P= -H =1 « i? 2 /2, R 2 = R/2A. (26) 
The correlation function at small q 2 limit, qR <C 1, is then: 

C(p, q) = 1 + ^4^4 - 1 + -P (~ ,"' (p) , ^ q 2 ) • (27) 



(a£a Pl ) (a£ 2 a P2 ) \ n g (p)+n c (p 



At large <? limit, qR 3> 1, it is 



2 
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C(p,g) = l+ — exp(-i^). (28) 



The effective interferometry radius squared, corresponding to C(q e ff) = 1 + 1/e, is 

^e// = <Utt = R ll t 1 + 2 ln K(p)/%(p))] • (29) 

Finally, we can see that when the density in phase space increases the interferometry radius of the gas component 
reduces to R/V% at most (see (26)). At very large densities, exceeding essentially the critical one, the Bose condensate 
determines the correlation function behavior at small p. Then the interferometry radius is reduced additionally as 
compared with its behavior for pure Bose-gas. At n/n cr — > oo, the interferometry radius in inclusive measurements 
tends to zero whatever large is the geometric size of the system. Note that intercept of the inclusive correlation function 
is equal to 2 and does not depend on n/n cr . This reflects the chaotic nature of the thermal multiboson source. We 
are not discussing here the very complicated problem of spontaneous symmetry breaking of the condensate in finite 
systems which could take place when there is the interaction destroying the degeneration in the system. 
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